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Evidence for color fluctuations in the nucleon in high— energy scattering 
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We study quantum fluctuations of the nucleon's parton densities by combining QCD factorization 
for hard processes with the notion of cross section fluctuations in soft diffraction. The fluctuations 
of the small-x gluon density are related to the ratio of inelastic and elastic vector meson production 
in ep scattering. A simple dynamical model explains the HERA data and predicts the x— and 
Q 2 -dependence of the ratio. In pp/pp scattering, fluctuations enhance multiple hard processes (but 
cannot explain the Tevatron CDF data), and reduce gap survival in central exclusive diffraction. 
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Hard processes in high-energy ep and pp/pp scattering 
probe nucleon structure at a resolution scale where it can 
be described in terms of the quark and gluon degrees of 
freedom of QCD. Essential in the analysis of such pro- 
cesses is the method of factorization, by which the am- 
plitude or cross section is separated into a short-distance 
quark/gluon subprocess, calculable in perturbativc QCD, 
and the distributions of the partons in the initial and final 
hadrons, governed by long-distance, non-perturbative in- 
teractions. Inclusive scattering experiments probe the 
longitudinal momentum densities of the partons. Mea- 
surements of hard exclusive processes in ep scattering re- 
veal information also about their spatial distribution in 
the transverse plane (generalized parton distributions, or 
GPDs). Such experiments can eventually provide us with 
a full 3-dimensional image of the quark/gluon single- 
particle structure of the nucleon. 

From the perspective of many-body physics the parton 
densities represent average characteristics of the nucleon, 
reflecting the quantum-mechanical average over config- 
urations in the nucleon wave function of different size, 
number of particles, etc. Equally fundamental are the 
fluctuations around the average value, which generally 
provide information about the nature of the interaction 
in the system. In the case of the quark and gluon densi- 
ties the fluctuations are related to variations of the size 
and intensity of the long-wavelength color fields inside 
hadrons — information crucial for understanding hadron 
structure in terms of non-perturbative QCD. An inter- 
esting question is which experimental observables could 
reveal such "color fluctuations" inside hadrons. 

To see how this problem might be approached, it is 
worthwhile to recall some facts about soft diffractive 
hadron-hadron scattering at high energies and low mo- 
mentum transfer, t. In such processes the hadrons can 
be represented as a superposition of states which diag- 
onalize the T-matrix and experience only individual at- 
tenuation ("diffractive eigenstates"). Elastic diffraction 
then results from uniform attenuation of these compo- 
nents, while inelastic diffraction is caused by different 



attenuation, which destroys the coherence of the original 
hadronic wave function. This idea was originally formu- 
lated using the language of the eikonal approximation 
appropriate for non-relativistic scattering P, Q ■ Later a 
more general formulation was developed, consistent with 
the fundamental principles of relativistic quantum field 
theory (causality, energy-momentum conservation). It 
uses the concept of the cross section distribution, P(cr), 
which describes the probability for the hadrons to scat- 
ter in a configuration with given cross section [2], |3| . It is 
normalized as / da P(a) = 1, and its variance is related 
to the ratio of inelastic and elastic soft diffractive cross 
sections at t = Q, 
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where the brackets denote the average with the distri- 
bution P(a). In particular, this formulation allows one 
to incorporate color transparency (the vanishing of the 
interaction of small-size configurations) in the a — > 
behavior of P(a), and the approach to the black-disk 
regime (unitarity limit) at high energies; both are funda- 
mental predictions of QCD. Detailed phenomenological 
studies of various diffractive phenomena involving pro- 
ton and nuclear targets have shown the usefulness of this 
approach and determined the properties of P(a) 0, 0j- 
At pp energies y/s ~ 20GeV the variance is u) a ~ 0.25, 
indicating sizable fluctuations of the interacting configu- 
rations in soft processes. 

In this Letter we propose to study color fluctuations in 
hadrons by combining QCD factorization for hard pro- 
cesses with the notion of cross section fluctuations in 
soft diffractive processes. We introduce the concept of 
a configuration-dependent parton density and follow its 
implications for various types of high-energy scattering 
experiments with hard processes. Our investigation pro- 
ceeds in three stages. First, we relate the fluctuations of 
the gluon density to the ratio of inelastic to elastic hard 
diffraction in ep scattering (HERA, future Electron-Ion 
Collider, or EIC). Second, we propose a simple model of 
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color fluctuations in the nucleon to illustrate and quantify 
our results. Third, we discuss the implications of color 
fluctuations for pp/pp collisions with multiple hard pro- 
cesses, and for rapidity gap survival in exclusive diffrac- 
tive scattering (Tevatron, LHC). The effects described 
here are not included in present Monte-Carlo generators 
for pp collisions. A more detailed account of our studies 
will be given elsewhere @. 

Consider diffractive production of vector mesons in ep 
scattering at Q 2 > few GeV 2 , t£ + p — ► V + X, where 
the proton may remain intact or dissociate into hadronic 
states X. The initial proton state can be expanded in a 
set of partonic states characterized by the number of par- 
tons and their transverse positions, summarily labeled as 
\ n ) : \p) — J2 n a n\n). Each configuration n has a definite 
gluon density G{x, Q 2 \n), given by the expectation value 
of the twist-2 gluon operator in the state \n), and the 
overall gluon density in the proton is 

G(x,Q 2 ) = J2 n K\ 2 G(x,Q 2 \n) = <G). (2) 

Because the partonic states appear "frozen" on the typi- 
cal timescale of the hard scattering process, one can use 
QCD factorization to calculate the amplitude for vector 
meson production configuration by configuration. It is 
(up to small calculable corrections) proportional to the 
gluon density in that configuration [7[ • An essential point 
is now that in the leading-twist approximation the hard 
scattering process attaches to a single parton, and, more- 
over, does not transfer momentum to that parton; it thus 
does not change the partonic state \n). Making use of the 
completeness of partonic states, we find that the elastic 
(X = p) and total diffractive (X arbitrary) cross sections 
are proportional to 



(da cl /dt) t=0 ex [EnK\ 2 G(x,Q 2 \n)] = (G) 2 , (3) 
(da diS /dt) t=0 oc £> n | 2 [G(x,Q 2 \n)] 2 = (G 2 ). (4) 



For the cross section of inelastic diffraction <7i„ c i = Odiff - 
<r e i we thus obtain 
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This model-independent relation allows one to infer the 
fluctuations of the gluon density from the observable ra- 
tio of inelastic and elastic diffractive vector meson pro- 
duction. It can be easily generalized to a large variety of 
hard processes such as 7^ + T — > 2n (two jets) + T, or T 
production in ultraperipheral pp collisions at LHC [8( . 

Generally, we expect w g to be a weak function of Q 2 
at fixed x (approximate scaling), as the gluon density 
depends only logarithmically on Q 2 . The x-dependence 
of uj g at fixed Q 2 is difficult to infer from first principles; it 
depends on the "color flow" in the nucleon wave function, 
i.e., how the small-x parton densities change with the 
configuration of the large-x constituents. 



0.3 



0.2 



0.1 





1 1 

Q 2 [GeV 2 ] = 3 




10 — 




100 — 


■ 


1000 ■ 







10- 



10" 



10" 



10" 



FIG. 1: The variance of fluctuations of the proton's gluon 
density, lj s , as a function of x for several values of Q 2 , as 
obtained from the scaling model, Eq. (|6]), and a phenomeno- 
logical parametrization of the gluon density. 



To estimate the variance of the gluon fluctuations at 
small x, and to study their implications in other hard 
scattering processes, we propose here a simple model 
based on two assumptions: (a) At moderate energies 
(y/s ~ 20 GeV) the hadronic cross section of a configura- 
tion is proportional to the transverse area occupied by the 
color charges in that configuration, a oc -R 2 onfig ; (b) the 
normalization scale of the parton density changes propor- 
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tionally to the size of the configuration, fi cx R : 
cr _1 . Assumption (b) is similar to the "nucleon swelling" 
model of the EMC effect |9| and implies a simple scaling 
relation for the cr-dependent gluon density: 



G{x,Q 2 \a) = G(x^Q 2 ), 

£(Q 2 ) = (*/(a)) a * (Q2 ° )/a ° (Q2) , 



(6) 



where Q\ ~ 1 GeV 2 . Assumption (a) then allows us to 
compute the configuration average using the phenomeno- 
logical cross section distribution found in Ref. Fig- 
ure [1] shows the result for the variance of the gluon den- 
sity in this model. At small x and low Q 2 it reaches 
values comparable to the variance of soft cross section 
fluctuations, u>. Note that our evolution-based model of 
gluon fluctuations applies primarily to small x (<C 0.1); 
at larger x non-perturbative correlations not included 
here may become important. Present experimental data 
on the cross section ratio in Eq. ([5]) are very limited. 
The value uj g ~ 0.15 - 0.2 for Q 2 = 3 GeV 2 and x ~ 
10~ 4 - 10~ 3 obtained in our model is consistent with the 
HERA data on vector meson production, where the ef- 
fective scale is Q 2 S ^2 — 4 GeV 2 . The data also indicate 
weak dependence of the ratio on Q 2 and the vector me- 
son mass; however, the limited Q 2 range and the lack of 
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dedicated studies do not allow us to test our model pre- 
dictions in more detail. Future measurements with LHC 
and EIC could significantly improve the situation. 

Correlations between fluctuations of the parton densi- 
ties and the soft-interaction strength have numerous po- 
tential implications for high-energy pp/pp collisions with 
hard processes. One example is the relative probability 
of double binary parton-parton collisions (see Fig. [2^), 
defined as the ratio 
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where O12 etc. are the variables characterizing the ob- 
served dijets (or photons), and f(x\) and f(x\,x 3 ) etc. 
are the single and double parton densities, respectively 
(we suppress the dependence on the scale) . The effective 
cross section <7 e ff is a measure of transverse correlation 
between partons. The FNAL CDF experiment measured 
a value of cr ff = 14. 5± 1.7^3 mb in 3-jet + 7 events [Icj . 
a factor of ~ 3 smaller than the naive estimate based on 
a uniform distribution of partons with the e.m. radius 
of the nucleon. Proper account of the transverse geom- 
etry, using gluon GPDs based on J/ip photoproduction 
data, reduces the discrepancy to a factor of ~ 2 11[ . An 
interesting question is whether fluctuations of the sizes 
of the interacting configurations and the gluon density 
could explain the remaining discrepancy. 

Elementary geometric arguments show that in the 
mean-field approximation (i.e., no fluctuations) 

a J (mean field) = Jd 2 b P 12 (b) P 3A (b), (8) 

where b = \b\ is the pp impact parameter, and P\ 2 {b) 
describes the ^-dependence of the probability to have 
two partons with x\ and x 2 collide in the pp collision, 
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= Jd 2 Pl Jd 2 p 2 S^(b 
x F(x 1 ,p 1 ) F(x 2 ,p 2 ), 
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FIG. 2: (a) Double hard scattering in high-energy pp- 
collisions. (b) Schematic illustration of fluctuations. In the 
model used here [cf. Eq. JBJ], configurations with larger size 
have larger gluon density. 



with J d 2 bPi 2 (b) — 1 (a similar definition applies to 
P34). Here F(x,p) is the normalized transverse spatial 
distributions of the partons obtained from the GPD. To 
estimate the effect of fluctuations, we now apply the 
same reasoning to the instantaneous configurations in 
the colliding protons, and let both the transverse size 
and the gluon density in these configurations fluctuate 
as described by our model [cf. Eq. (J6])]. We assume 
that the transverse size of the configurations is pro- 
portional to a, i.e., the spatial distribution of partons 
is F(x,p\a) = \~ 1 F(x, \~ 1 / 2 p), with A = a /(a) (see 
Fig. Ht>) ■ With these assumptions we obtain 

a~g (fluct) = (( (/)(xi\ai) 4>(x 2 \a 2 ) (f>(x 3 \ai) (f>(x 4 \a 2 ) 
x fd 2 b P 12 (b\ai,a 2 ) PuibWuaz) )) 12 , 
<Kzi|<7i) = /(zi|ai)/</(a;i|<7i)}i etc., (10) 

where the double brackets denote the average over the 
cr-distributions for both protons. 

Because in multijet production at the Tevatron the 
typical x-values of the partons are large (~ 0.1), fluc- 
tuations of the parton densities are much smaller that 
in vector meson production at HERA (see Fig. [IJ. The 
dominant effect in Eq. (fTOf thus comes from fluctuations 
of the sizes of the interacting configurations. Keeping 
only the latter, we obtain a simple analytic result by ex- 
panding in leading order in the variance uj a : 

cr c ff (fluct) = (1 — CT ff (mean field) ; (11) 

numerical studies show that higher-order corrections are 
negligible in practice. One sees that size fluctuations 
indeed reduce a e g, because of the disproportionate en- 
hancement of multiple hard processes in small-size con- 
figurations. However, the reduction in our model is found 
to be only of the order 10 — 15%, which cannot account 
for the discrepancy with the CDF value. This indicates 
that other dynamical mechanisms must be responsible 
for the enhancement of multi-parton collisions, e.g. local 
transverse correlations between partons as suggested by 
a "constituent quark" picture of the nucleon [12J . 

We note that much stronger reduction of a e s was ob- 
tained [H] in a simplified multi-channel eikonal model, 
where all hadronic diffractive states are represented by 
a single channel [3]. In such a model, to fit the avail- 
able data on elastic, inelastic, single and double diffrac- 
tive cross sections, one needs to enhance the strength of 
the Pomeron coupling between diffractive eigenstates of 
smaller radii. In the present model this would correspond 
to an increase, rather than a decrease, of the parton den- 
sity for small a, which is the origin of the different results. 

Another class of processes affected by fluctuations is 
central exclusive diffraction, pp — > p + H + p (H = dijet, 
heavy quarkonium, Higgs boson), widely discussed as a 
possible channel for the Higgs boson search at the LHC. 
Of interest here is the rapidity gap survival (RGS) proba- 
bility, S 2 , measuring the suppression of the cross section 
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FIG. 3: Fluctuations in central exclusive diffraction pp — > 
p + H + p. One of the protons is treated as a superposition of 
partonic states |n), which undergo diffractive scattering with 
the other proton (small fluctuation approximation). 

resulting from the requirement of no inelastic soft spec- 
tator interactions. This quantity has extensively been 
studied using models based on eikonalized Pomeron ex- 
change [HI]. A recent analysis in a partonic approach 
with account of the transverse geometry (GPDs) found 
that the suppression of the diffractive cross section results 
mainly from the elimination of pp collisions at small im- 
pact parameters, in which there is a high probability of 
inelastic interactions (l6j . This result was obtained in the 
mean-field approximation, in which the GPDs and the 
soft-interaction strength are taken at their average val- 
ues, with no correlations between them. (The approach 
of Ref. [IH gives results comparable to the mean-field 
approximation and also ignores such correlations.) 

To estimate the effects of correlations, we follow the 
approach described above and allow for configuration de- 
pendence of both the gluon GPDs (i.e., the gluon den- 
sity and the radius of the transverse distribution) and 
the soft-interaction strength. The latter we model phe- 
nomenologically, using data on the energy dependence of 
(a) and oj a in soft hadron-hadron scattering [6[ . At small 
x the dominant effect comes from the correlation of fluc- 
tuations of the gluon density with the soft-interaction 
strength. Treating the fluctuations as a small correction 
to the mean-field result we have 

S 2 (fluct) = (l + 4e) S 2 (mean field), (12) 

where e is the correction obtained if one of the protons 
fluctuates (see Fig. [3]); the factor 4 counts the number 
of protons in the diffractive amplitude and its complex 
conjugate. We find e w —0.07 (—0.04) for production of 
a system with mass Mh = 10 (100) GeV at -y/s = 2 TeV 
(Tevatron), resulting in a reduction of the RGS proba- 
bility by a factor ~ 0.7 (0.85). The sign reflects the fact 
that smaller configurations with higher survival probabil- 
ity have a lower density of small-a; gluons in our model 
(see Fig. [2}d). We note that at higher energies (LHC) 
the onset of the black-disk regime in hard interactions 
causes another, more substantial reduction of the RGS 
probability relative to the mean-field result fill ]. 

Fluctuations also affect the final-state transverse mo- 
mentum dependence of the cross section for central ex- 



clusive diffraction, because they change the size of the 
dominant interacting configurations. Our findings sug- 
gest that the px distribution at RHIC and Tevatron en- 
ergies is narrower than given by the mean-field approx- 
imation — a prediction which could be tested by future 
measurements of diffraction at RHIC. 

In sum, the study of quantum-mechanical fluctuations 
of the quark/gluon densities is the natural next step in 
the exploration of nucleon structure in QCD, following 
the mapping of the longitudinal momentum and trans- 
verse spatial distributions of partons. Detailed measure- 
ments of diffractive vector meson production (HERA, 
EIC) could significantly enhance our knowledge of fluc- 
tuations of the small-x gluon density. They also provide 
essential input for modeling the dynamics of high-energy 
pp collisions (RHIC, Tevatron, LHC), where fluctuations 
play an important role in multijet production and ra- 
pidity gap survival. How the concept of parton density 
fluctuations developed here is affected by the approach to 
the unitarity limit at high energies remains an interesting 
problem for further study. 
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